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ABSTRACT 


A high resolution 20-level, primitive equation model is initialized using sigma-0 
(O,) data acquired during a Coastal Transition Zone (CTZ) cruise from 6 to 12 July 1988 
near Point Arena, California for the purpose of diagnosing the three-dimensional 
dynamically balanced flow field for the region. The major feature in the region during 
the cruise period was a strong meandering jet which flowed equatorward and offshore, 
oriented in a northeast to southwest direction. A density diagnostic version of the model, 
in which the 0, field is held constant during the 72 hour integration time, is compared 
with a robust diagnostic version of the model, in which the 6, field is allowed to adjust 
dynamically within the constraints of a damping time scale. The density diagnostic 
version of the model produces vertical velocities which are an order of magnitude larger 
than vertical velocities calculated from bio-optical data collected in the same location and 
at the same time as the CTZ cruise. The robust diagnostic version produced vertical 
velocities within the range of vertical velocities calculated from the bic-optical data. The 
horizontal pattern of vertical velocity from the robust diagnostic version also compared 
well with patterns seen in trough/ridge systems in the Gulf Stream. In particular, there 
is horizontal convergence and downwelling all along the offshore directed part of the 
meandering jet. Vertical cross-section comparisons between model velocity data and 
ADCP velocity data collected during the cruise revealed differences which may be caused 
by factors neglected by the model (e.g. bottom topography, wind forcing, and heat fluxes). 
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I. INTRODUCTION 


A. BACKGROUND - THE CALIFORNIA CURRENT 

The California Current is the southeastward extension of the North Pacific Drift 
which flows across the Pacific Ocean from west to east at about 45° N (Tchernia, 1980). 
This current flows toward the equator, and is considered a cold, easter boundary current. 
It. along with three other currents, comprises the California Current System. The other 
three currents are: the Califomia Undercurrent, the Davidson Current and the Southem 
California Countercurrent (Hickey, 1979). These currents are all poleward flowing, the 
Davidson Current and the Southern California Countercurrent being surface currents, and 
the Califomia Undercurrent flowing beneath the surface. The Davidson Current is found 
north of Point Conception during the fall and winter. The Southem California 
Countercurrent exists to the south of Point Conception in the Califomia Bight. 

The Califomia Current System impacts the west coast states climatically and 
economically. Upwelling associated with it brings colder, nutrient rich water to the 
surface, supporting marine life there (Brink, 1983; Huyer, 1983). The colder water at the 
surface also cools the relatively warm air flowing over it, producing mild weather 
conditions and fog along the coast (Foster, 1989). 

Cold filaments, which occur in the California Current System from Vancouver 
Island to Point Conception. are of particular interest. These filaments consist of cold. 


salty inshore water which has been carried as far as 300-400 km offshore from the coast 











(Strub and the CTZ Group, 1990). Infrared and visible color satellite images reveal these 
filaments in the California Current System as well as in other eastern boundary currents 
(Traganza et al, 1983; Abbott and Zion, 1985; Flament et al., 1985; Kelly, 1985). 

The Coastal Transition Zone (CTZ) experiment examines some of these cold 
filaments in detail. Its purpose is to study the physical, chemical and biological structure 
of the filaments which occusted off the northem Califomia coast during the spring and 
summer of 1987 and 1988 (Strub and the CTZ Group, 1990). The goals of the CTZ 
experiment are: understand the interaction processes between the cold filament water and 
the open ocean water; study the effect the filaments have on the biology and metereology 
of the coastal zone; develop maps of the current, hydrographic, nutrient, and particle 
fields; and use numerical models to investigate what driving forces, dynamics, topography 
and stratification is necessary to duplicate the observed characterisitcs of these filaments 


(CTZ Group, 1988). 


B. OBJECTIVE 

The objective of this thesis research is to diagnose the velocity fields obtained from 
a 20-level robust primitive equation model using sigma-8 (o,) data from a cruise made 
during the 1988 CTZ field experiment. The cruise domain is from approximately 37°10’ 


to 39°21’N and 123°50’ to 126°40’W, and covers about one meander of the California 


Current. 








The same domain and model, lacking the robust feature, were used previously to 


diagnose the velocity fields from temperature data collected during the same experiment 


(de Jesus, 1990). 





II. DATA ACQUISITION AND PROCESSING 


A. DATA ACQUISITION 

Data from a Naval Postgraduate School cruise made from 6 to 12 July 1988 on the 
R/V POINT SUR that was part of the CTZ experiment are used here. This cruise is 
subsequently referred to as CTZRI. 

Acoustic doppler current profiler (ADCP) data and conductivity, temperature and 
depth (CTD) data were collected during the cruise. CTD measurements were made with 
an NBIS CTD system. Temperature measurement accuracy and salinity measurement 
accuracy were about +0.01 C and £0.003 psu respectively. Data were collected to a depth 
of 500 m. The ADCP made by RD Instruments operated continuously during the cruise. 
The transducer had a nominal range of 350 m and operated at 150 kHz with a 4 m pulse 
length. Retuming echoes were sampled selectively in 4 m bins to produce a vertically- 
averaged, ship relative velocity. Relative velocity profiles were made every three 
minutes. Reliable data were collected to depths as shallow as 25 m (Huyer et al., 1990, 
de Jesus, 1990). 

The horizontal sampling grid spacing is nominally 25 by 40 km (Stanton et al., 
1990). The actual station grid extends approximately 120 nm offshore and 110 nm 
alongshore, roughly parallel to the Califomia coast and southwest of Point Arena (Figures 
1 and 2). The alongshore sections are approximately 20 nm apart and have individual 


stations approximately 13 to 14 nm apart (de Jesus, 1990). Not all of the stations in the 





Figure 1. CTZ Station Grid for 1988: (from Huyer et al., 1990). 


CTZ grid are included in the model domain to eliminate coastal effects (Figure 2) (de 


Jesus, 1990). 


B. DATA PROCESSING 

The sigma-9 (6,) data were derived from the CTD data and objectively analyzed 
using the optimal interpolation method described, for example, by Carter and Robinson 
(1987). This method and some of the following procedures were similar to those used 
by de Jesus (1990) for temperature data, but with some differences as described below. 

The spatial mean 0, profile for the cruise and, the objectively analyzed 0, at each 
of the points in a grid with 4 km horizontal spacing was obtained with a vertical spacing 
corresponding to the upper 7 levels of the numerical model as given in Table lt. This 


mean 6, profile was calculated at meter depth intervals in the upper 500 m, corresponding 


ff Eureka 
sie Cape Nendocino 


ool yaw 130W 126W 126W 124W 122W 





Figure 2. Actual Topography in the Model Domain: The box shows the area of the 
model domain. Depths are in meters. 


to the depths at which observed data were collected, and then averaged in the vertical to 
every 10 m. 

The mean 0, profile was extended to 4000 m by fitting an exponential curve to the 
mean profile in the upper 500 m using a centered difference scheme around 450 m. The 


6, at 4000 m was left to be determined by the exponential curve, increasing the degrees 








of freedom and allowing continuity of 0, and its vertical derivative at the interface 
between the observations and the extrapolation (i.e. at 450 m). 

This extended mean 6, profile was used to calculate the first three dynamical 
vertical modes 6,,(k), Oe,(k), and 04,(k) following the procedure described in Rienecker 
et al. (1987) and outlined below. 


The initial definition made was the perturbation 


Glick) = Oyliv,k) - 5, (A), {1} 


where i, j, and k are the gridpoint indices for the x, y, and z directions, respectively, 
Gy (i,j,k) is the perturbation part of 0, at the gridpoint, 6,(i,j,k) is the objectively analyzed 
G, at the gridpoint, and 6, (k) is the mean 6, at level k. 

To extend the perturbation below 500 m, 0,’(i,j,k) is assumed to be a linear 


combination of the first three dynamical vertical modes, such that 


Galiv.k) = 4,(ij)0,,(k) + a,(if)O,(k) + 4,(i 0,4), {2} 


where a,(i,j) is the weighting given to 6,,(k), a,(i,j) is the weighting given to 6,,(k) and 
a,(i,j) is the weighting given to 6,,(k). Solutions for a,, a,, and a, at each horizontal grid 
point i,j need to be calculated. 

The procedure for determining a,, a,, and a, followed by de Jesus (1990) involved 
solving a system of four equations in three unknowns (a,, a,, and a,) using a linear least 


squares method. This procedure. referred to as the "mean conserving” formulation, was 





also used here. In addition, two new procedures were used which gave different results 
for the vertical profiles below 500 meters. They were considered possibilities for better 
representing the physical condition of the ocean in the cruise domain. 

The "mean conserving" formulation of the linear least squares method is outlined 
briefly as are the two new, possibly more physical, formulations. For the sake of 
comparison, vertical profiles of 6,’ are calculated at representative (i,j) points using each 
of the three methods. Arguments are made conceming which of the three methods are 
considered the best for representing the characteristics of the ocean. Finally, an exact 
method of solving for a,, a,, and a, using four equations and four unknowns is compared 
with the linear least squares method, and arguments are made for using the linear least 
squares method instead of the exact method. These issues are discussed further below. 

Following the "mean conserving" formulation of the linear least squares method, 
four equations are derived from equation {2}: 

<Oa(ij,k)> =a,<0,,(k)> +,<0,,(k)> +4,<0,,(k)>, {3} 
<O,,(K)06(i,k)> =4,<0,,(K)O, (k)> +4,<0,,(k)0,,(k)> +4,<0,,(K)0,,(K)>, 4} 
<0, (10 a(iy,K)> =4,<0,,()0,, (> +4,<0,,(b)0,,()> +4,<0,,(K)0,,(K)>, 5} 


<64,(K)Oo(i,k)> =4,<0,,(K)G,,(k)> +4,<0,,(k)0,,(k)> +4,<0,,(k)0,,(k)>. {6} 


The bracketed (< >) terms are a vertical average, defined as 


7 7 
<a>= PY QAnH/yY A xh. (7) 
1 1 





where Q is any expression inside the brackets, Az(k) is the vertical thickness between 
adjacent k levels in the model, and the limits 1 and 7 with the summation sign correspond 
to levels | and 7 in the model(Table 1). Therefore, <Q> is the vertical average of Q(k) 
over the top seven levels of the model, which are encompassed by the analyzed 
observations (i.e. the upper 500 meters). Equation (3} is the statement of “conservation 
of the mean” in this formulation. 

This system of four equations in three unknowns (a,, a,, a) is solved by the linear 
least squares method. This means that none of the four equations are satisfied exactly, 
but that there is a residual error in each one. After a,, a,, and a, are found for each i,j 
point in the domain of the model, equation {2} is used to find 6, (i,j,k) for k = 8 tok = 
20 (i.e. from 500 meters to 4000 meters) (Table 1). From Equation {1}, 04(i,j,k) is found 
for all 20 levels in the model. 

The first of the two new formulations alluded to above represents an energy-related 
approach to arriving at solutions for a,, a,, and a;. It involves the replacement of equation 


(3) by 


<0 ij IO (i K)> =4,<O4(i KO, (K)> +4,<04(i Ky, (b)> +4, <04(i,K)G,,(4)>- {8} 


The 0,” on the left hand side of the equation is an energy term. Since this equation puts 
a constraint on the energy of the 0, in the upper ocean where the observations occur, this 


is called the "energy conserving” formulation of the problem. 


The second of the new formulations puts a constraint on the continuity of 0, at the 
interface between the observations and the region of extrapolation below. In this case. 


the interface occurs at level 7 in the model. Therefore, equation {3} is replaced by 


O9(iyf,7) =a,5,,(7) +4,0,,(7) +4,0,,(7). {9} 


Because of the constraint at the interface, this is called the "interface matching” 
formulation. Since both of the new formulations are also solved using a linear least 
squares method, there is a residual error in each equation as there is with the first 
formulation. 

When comparing the results of vertical profiles using these three formulations, 
differences among them can readily be seen. In the CTZR1 cruise domain, two points 
where 6, is large in the upper levels of the ocean are used to illustrate these differences 
(Figures 4 and 5). 

The (i,j) gridpoint (27,36), which is approximately map coordinate (x,y) point 
(110,144) in Figure 3, is located in a region of low density water, with 6, < 24.4 kg/m’ 
at 65 m. Figure 4 shows the vertical profile of 0, at this point for the three different 
formulations of the linear least squares solution to a,, a,, and a, in equation (2). Clearly, 
the observed profiles in the upper 500 meters are the same in each case. Also, below 
about 2000 meters the profiles are the same. In the deep ocean, the water characteristics 
are such that the dynamical modes are all small in amplitude so that changing the 
formulation of the solution does not affect the profile. Between these two levels, the 


three formulations produce three different profiles. The 6, profile for the "interface 
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Figure 3. Plot of o, at 65 Meters Depth(Model Level 3), Cruise CTZR1, 06-12 
July 1988: The two points of interest are marked. 


matching” formulation has the smallest amplitude of all. The constraint on having 
continuity at the interface between the observations and the extrapolation is the cause of 


this reduced amplitude. The "energy conserving” formulation is also smaller in ainplitude 
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than the "mean conserving" formulation of the solution. Comparing the profiles in Figure 


4 with the three vertical dynamical modes shown in Figure 6 suggests that mode 2 and 
mode 3 are the main contributors to the structure of the vertical profile of 0,’ at the level 
between 500 and 2000 meters. 

The (i,j) gridpoint (21,8), which roughly corresponds to map coordinate (x,y) point 
(85,32), is in an area of high density water at 65 m as shown in Figure 3. The vertical 
profiles of the three formulations at this point are not as distinctly different as those at 
point (110,144) as seen when comparing Figure 5 with Figure 4. At point (85,32), the 
“interface matching” formulation and the “energy conserving” formulation have a smaller 
amplitude than the “mean conserving" formulation in the level between 500 meters and 
2000 meters. This illustrates the differences among the points in the cruise domain 
resulting from the 3 different extrapolation formulations. On the other hand, the 
similarities of the profiles above 500 meters and below 2000 meters can be explained 
using the arguments made previously for point (110,144). 

Of the three formulations described above, which is the "best" to use for input into 
the model? The "mean conserving” formulation suffers because it places no constraint 
on the amplitude of o,. The best example of this is a sine wave. The mean of a sine 
wave is the same no matter what its amplitude. The “energy conserving” formulation 
may alias high wave number energy observed in the upper 500 m (which actually belong 
to higher dynamical modes) onto the lower wave number modes, which is due to the 
energy dependence caused by the squared term on the left-hand side of equation (8}. 


This would tend to increase the amplitude of 0,’ compared to a method that did not alias 


14 











‘sishyeue sansefgo oy) 
WO paye[nofes °o Jo apyoid Jeonsea ues TeNeds popuayxe oy) WIJ 


pandwod :FWZLD 105 °O jo sapoyy yeormeudg yeonssA sosy yy psy “9 aanBry 


ag EE AS, 
MODE 3 


cH/94 NI @ UWIIS 


1500.0 2000.0 2500.0 3000.0 3500.0 4000.0 


DEPTH IN METERS 


Oo 
9 
O 
~ 


500.0 


0.0 














the high wave number variability into the solution. Therefore, the “interface matching” 
formulation seems to be the "best" to use for the model. As stated above, it places a 
constraint on the 0, at the interface between the observed data and the extrapolated data, 
and does not alias high wave number variability onto the three low wave number modes. 
The problems associated with the other two methods are, therefore, eliminated. However, 
before a definite conclusion can be made about the relative ability of these formulations 
to depict sub-500 m variability, they should be tested with full water column observed 
data. Because of the arguments given above, the "interface matching" formulation is used 
in this study. Maps of a,, a,, and a, using the "interface matching” formulation are shown 
in Figures 7, 8, and 9. These modal amplitudes have large spatial coherence, and their 
patterns closely correspond to the 6, field in Figure 3. 

The problem of solving for a,, a,, and a, in equation {2} could have been done 
exactly with four equations and four unknowns using variational techniques found in 
Hildebrand (1965). In Hildebrand (1965; see pp. 119-121), it is shown that the problem 
of determining a,, a, and a, that minimizes the root mean square difference between the 
observed 6,’ and the 6, given by (2} in the upper 500 m subject to the constraint (9}, 
can be expressed as an unconstrained minimization problem which leads to the following 


four equations 


16 





PEO bh ib. Vilage» 
CRUISE, CIZKI 5: O6=12° JULY 


ens 

wou? 

maw eh 
wewanay 
ne oe 


em 


way 


“Var Z eens = 


Oo 
N 
N 
Oo 
Oo 
N 
oO 
oO- 
-_ 
oO 
w 
—_— 
fen} 
ba 
od 
Oo 
Nn 
-_ 
fen) 
oO 
— 
oO 
oO 
oO 
wo 
Oo 
~~ 
Oo 
N 
oO 


0 20 40 60 80 100 120 140 160 180 200 22U 240 260 280 300 
X 





Figure 7. Plot of A, Values, Cruise CTZR1: Obtained from "interface matching” 
formulation described in text. 
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Figure 8. Plot of A, Values, Cruise CTZR1: Obtained from "interface matching" 
formulation described in text. 
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Figure 9. Plot of A, Values, Cruise CTZR1: Obtained from "interface matching’ 
formulation described in text. 
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Delii,7) =4,09,(7) +4,09,(7) +2,04,(7)s {10} 


<0y, (A) Oe(i,.k)> =4,<6,, (0, (A)> +4, <0, ,(k)0,,(k)> +4,<0,,(K)0,,(k)> -20,) {11} 
<q, (KG ali, k) =A, <0y,(K)Dy,(K)> +4, <q, (KIO, (K)> +4,<g,(k)Og5(K)> - F04(7), {12} 
<0,,(k)04(i,k)> =0,<0,,(4)0,,(K)> +4,<0,,(1)0,,(K)> +4,<0,,(4)09,(K)> - *0,17) {13} 


where A is the (unknown) Lagrange multiplier. Equation {10} is the same as equation 
(9}. Equations (11), (12}, and (13} are the same as equations {4}, (5}, and {6}, 
respectively, with the exception of the final term on the right hand side of the equations. 
This final term, the term which includes the Lagrange multiplier, is a residual term. 
Notice that this exact method allows the constraint, in this case {10}, to be satisfied 
exactly. All the residual error is put into the other 3 equations. As pointed out above, 
the least squares method distributes the error among all 4 equations, none of which are 
satisfied exactly. These residuals represent the extent to which 6, as represented by the 
first three vertical dynamical modes does not match the actual 6, in the region of the 
ocean where actual observations have been made, in this case the upper 500 m. The 
residual term, therefore, represents a measure of the high wave number variability above 
500 m that simply cannot be represented by a sum of the first 3 vertical dynamical 
modes. Since 0,’ may have high wave number variability at all levels at and above 500 


m, it does not seem physically realistic to force an exact contraint at level 7 of the model 


as in equation {10}. 








Thus, the least squares method, with 4 equations and 3 unknowns, was chosen for 
this study. In this least squares method, none of the equations was satisfied exactly, but 
the total residual, spread over all 4 equations, may be smaller than it would be using the 


exact method. This is also a subject worthy of much further study. 


C. THE MODEL 

The numerical model is a 20-level, primitive equation model of a baroclinic ocean 
on an f-plane, with hydrostatic, Boussinesq, and rigid lid approximations, and a flat 
bottom (Haney, 1985). The space-staggered B-scheme is used for horizontal finite 
differencing, and a sigma coordinate system defines the vertical dimension (Arakawa and 
Lamb, 1977). The boundary conditions at the sea surface consist of no wind stress and 
no buoyancy flux. The lateral boundaries of the model domain are treated as open 
boundaries. A sponge layer of 5 gridpoints based on the discussion of Lorenzetti and 
Wang (1986) is used at the lateral boundaries around the entire model domain to 
eliminate the effects of gravity wave noise there. The coefficients of the sponge at the 
gridpoints, starting from the boundary inward, are .20, .13, .07, .03, and .008, 
respectively. The sponge affects only the velocity fields in the model. The total depth 
of the ocean is 4000 meters. Model levels are listed in Table 1. The u and v components 
of velocity, and 6, are calculated at the model levels, but the w component of velocity 
is calculated at the midpoints between the model levels as a vertically staggered 


arrangement. The values of constants used in the model are listed in Table 2. 
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The CTZR1! crise domain included 76 objectively analyzed gridpoints in the east- 
west direction by 60 gridpoints in the north-south direction. The model domain included 
only 64 of the gridpoints in the east-west direction by 60 gridpoints in the north-south 
direction, with the 12 eliminated gridpoints in the east-west direction being in the 
eastern part of the cruise domain along the California coast. This makes the model 
domain more uniform, and reduces any shallow water effects which may occur. 

Since the model grid points are much closer together than the observed data points, 
the objectively analyzed 0, field was somewhat noisy on the fine 4 km model grid. Since 


this noise is not considered to be physical, the following horizontal smoother was used: 


of = 55 (4m +8 +0), {14} 


where Ogy = Og (i,j,k), the perturbation 0, at each domain gridpoint, 

Og” = Oe (i+1, j,k) + Oe (i-1,j,k) + Oe (i,j+1,.k) + 0, '(i,j-1,k), and 

Oe" = Oo (i+1,j+1,k) + Og (i-1 j+1,k) + Og (i-1,j-1,k) + O,(i+1,j-1,k). This smoother was 
applied twice. Figure 3 is an example of a smoothed o, field. The corresponding 
unsmoothed field is shown in Figure 10. 

Two methods were used to diagnose the currents from a given 0, field. In both 
methods, the 0, field and the geostrophic currents (having zero vertical average over the 
full model depth) computed from the hydrostatic pressure field were used to initialize the 
20 level model. The model was then integrated forward in time (for several days) until 
the currents and density (0,) fields came into a new "balance." In one method. the 


det.sity diagnostic method. the 06, field is held constant throughout the model run. thus 


22 


Table 1. K-LEVELS AND CORRESPONDING DEPTHS IN METERS USED IN 
MODEL. 


2(k) (meters) 2k) (meters 





requiring the velocity fields to adjust to it. In the second method, a robust diagnostic was 
added, in which the 6, field was allowed to adjust, but was damped back to the initial 
features on a time scale which was short compared to the advective time scale, but long 
compared to the geostrophic adjustment time scale (Sarmiento and Bryan, 1982). This 
was done to hold the larger scale features constant, therefore, maintaining the diagnostic 
feature of the model, but allowing the smaller scale features to adjust. Much of the small 
scale noise arises from the multi-level initialization from the objectively analyzed, sparse 
observations. The scaling arguments used to justify the choice of the damping time scale 


are presented in the results section below. 
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Ill. MODEL RESULTS AND DISCUSSION 


The primitive equation model was initialized with the o, fields at the levels 
described above, and with geostrophically balanced currents having a zero vertical 
average. The model was then integrated forward in time for 72 hours. In the density 
diagnostic version of the model, the 0, field was held fixed. During this integration, the 
currents underwent adjustments to a new state of balance in which non-linearity and eddy 
diffusion effects were not negligible because the currents were no longer in exact 
geostrophic balance. Time series of the velocity components are used to show the process 
by which this balance occurs. The resulting diagnostic currents are discussed below. 

In addition, the robust diagnostic version of the model, as briefly described above, 
was integrated forward in time to 72 hours. The effects of the robust diagnostic version 
versus the density diagnostic version are discussed. Special attention is given to the 
vertical component of the velocity field and its relationship to vorticity advection for the 
robust diagnostic version of the model. Finally, a comparison is made between the 
currents as computed by the robust diagnostic version of the model and observed currents 
measured with an ADCP during CTZR1. 

To simplify the presentation and discussion of the results, only horizontal maps and 
time series at 65 m depth are used, which is a model level representing surface currents 
in the thermocline. Discussion of change in the vertical is limited to the comparison 


between ADCP results and model results. 
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Table 2. CONSTANTS USED IN MODEL 


fa feathers we i any 
C's a CR 
fax | ewridional gid wacng idem 
fy [seni gidpacng [een id 
[D___| total ocean depth 
LAT reference latitude 


biharmonic momentum dif- 4x 10" cm‘ sec" 

fusion coefficient 

biharmonic heat diffusion 4x 10" cm‘ sec! 
coefficient 















vertical eddy viscosity 
vertical eddy conductivity eee 


thermal expansion coeffi- 2.877 x 10* C' 
cient 


A. TIME SERIES OF VELOCITY FIELDS FOR DENSITY DIAGNOSTIC 
Was a new state of balance achieved by integrating the density diagnostic version 
of the model ahead 72 hours? Time series of the u, v, and w-components of velocity at 
(i,j) gridpoint (44,25) (ie. map coordinate point (172,96) in Figure 11) and at a depth of 
65 meters (90 meters for w because of the vertically staggered grid described in Chapter 
II) show the adjustments in the velocity at these locations (Figures 12 and 13). Note the 


large oscillations which occur during the early integration times in the model. These 
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Figure 10. Unsmoothed plot of 64. 


oscillations are damped as the new state of balance is approached at around 48 hours. 
The small oscillations continue to occur at 72 hours, especially in the u and v component 
velocities (Figure 13). The oscillations are remnants of the initial adjustments in the 


velocities, and could possibly continue for several more days. 
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This process can also be seen in the time series of the u, v, and w-components of 
velocity at (i,j) gridpoint (41,24), or (x,y) map coordinate point (160,92) in Figure 11. 
also at a depth of 65 meters (90 meters for w)(Figures 14 and 15). At this point, the 
damped oscillations beyond 48 hours in the u and v velocity components are almost 
sinusoidal. 

Examination of the time series of the component velocities at these two points show 
that an exact balance has not been reached in the density diagnostic version of the model 
at 72 hours. Small oscillations are still ocurring, although it appears that the currents are 
quite close to their final equilibrium values at 72 hours. The similar evolution of the 
velocity components with time at the two points indicates that the entire velocity field 
probably evolves in a similar way. Integrating the model beyond 72 hours to determine 
when the oscillations in the velocity components are fully damped was deemed 


unnecessary. 


B. SIGMA-@ AND VELOCITY FIELDS FOR DENSITY DIAGNOSTIC 

The ©, field at time 0 for a depth of 65 meters is shown in Figure 16, and the 
resulting hydrostsatic pressure field in Figure 11. By definition, these fields do not 
change with time in the density diagnostic model. 

The initial (geostrophic) u and v components of the velocity field at 65 meters are 
shown in Figures 17 and 18. The resulting vector velocity field is shown in Figure 19. 


Note the strong south to southwestward jet flowing from the northeastem comer to the 


southwestern comer of the model domain. Velocities in this jet reach as high as 80 
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Figure i1. Pressure Field of the Density Diagnostic at 65 Meters and 72 Hours: 
Points referred to are marked. Contour interval 2.5 cm. Dashed 


lines are negative values. 


cm/sec. Velocities as high as 75 cm/sec were observed at 50 meters by Huyer et al. 
(1990) during CTZR1. The velocities shown in the model domain at time 0 are simply . 


the geostrophic velocities with which the model is initialized. The initial w velocity 


component is zero by definition. 
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Figure 12. Time Series of W Velocity Component at Map Point (172,96) and 
Depth of 90 Meters for Density Diagnostic. 


After integrating the model ahead 72 hours, the velocities are no longer geostrophic. 
The u component velocity field has intensified, especially in the south central part of the 


model domain where the jet turns westward(Figure 20). The v component velocity field 
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Figure 13. Time Series of U and V Velocity Components at Map Point (172,96) 
and Depth of 65 Meters for Density Diagnostic. 


has also adjusted(Figure 21). This intensification and adjustment can also be seen in the 


time series shown in Figure 13. Note the increased u velocity shown at 72 hours as 


compared to 0 hours. Also note how the v velocity has changed. 


30 







DIME: SER TES “OF 


W VELOCITY (CM/S) 


0.06.0 12.0 18.024.0 30.0 36.0 42.0 48.054.060.066.072.0 
TIME CHRS) 


Figure 14. Time Series of W Velocity Component at Map Point (160,92) and 
Depth of 90 Meters for Density Diagnostic. 


The resultant horizontal velocity vectors for 72 hours are shown in Figure 22. A 
comparison of a few of the vectors shown here with the corresponding initial geostrophic 


velocities in Figure 19 reveals significant direction shifts as well as speed changes. The 
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Figure 15. Time Series of U and V Velocity Components at Map Point (160,92) 
and Depth of 65 Meters for Density Diagnostic. 


direction shift varies from point to point as illustrated by the time series in Figures 13 and 
15. From the u and v components of velocity in Figure 13, it can be inferred that the 


current has shifted westward by 72 hours. The direction of the current at the point in 
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Figure 16. Sigma-6 Field at Time 0 and Depth of 65 Meters: Contour interval .2 
kg/m? 





Figure 15 seems to have shifted to the east. 

The maximum velocities shown in the jet at 72 hours are over 100 cm/sec. This 
is higher than the maximum observed values noted above. 

The w component of velocity at 72 hours and 90 meters is shown in Figure 23. 
Maximum upwelling is around .8 cm/sec, and maximum downwelling is around 1.0 


cm/sec. At the points shown in the above time series, upwelling is occurring (Figures 12 
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Figure 17. U Component Velocity Field at Time 0 and Depth of 65 meters: 
Contour interval 20 cm/s. Dashed lines are negative values. 





and 14). At map point (172,96) upwelling velocity is around .25 cm/sec (Figure 12), and 
at map point (160,92) upwelling velocity is around .6 cm/sec (Figure 14), at 72 hours. 
These velocities convert to 200 to 700 meters per day which are much higher than the 6 
to 25 meters per day observed by Washbum et al. (1990) and Kadko et al. (1990) in the : 


CTZ region. 





34 





Cruise CTZR1, O6—12 July 1988 


VV hour 0, depth 65 m 





= 


. 
. 
e 
ere 
een 





weeene, 
~ 


. 





7 
‘. 
ny 










Pwene 
mewwe 


.. 
we esee, 










Pid wee 





. 
. 
- 


144 





Smee worn en 


128 


=r 
poneese . 
° 


342 


96 


Distonce (km) 





ete me wen, 








18 32 48 64 80 98 112 128 144 160 176 192 208 224 240 
Distance (km) 


0 
0 






Figure 18. V Component Velocity Field for Time 0 and Depth of 65 Meters: 
Contour interval 20 cm/s. Dashed lines are negative values. 


Therefore, the density diagnostic version of the model does not give horizontal or 
vertical velocities which agree closely with those observed. One possible reason for the 
rather poor results for the density diagnostic method is that the currents must adjust 
entirely to the given density field. However, in the real ocean there is a continual mutual 
adjustment between the currents and the density. Futhermore, on scales of motion that 


are sinaller than the intemal Rossby radius (about 30 kin in the CTZ region) the density 
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Figure 19. Resultant Horizontal Velocity Vectors at Time 0 and Depth of 65 
Meters: Vectors are plotted at every 4th gridpoint. 


field would actually adjust to the currents. Thus, the density diagnostic method is not 
really valid for scales of motion of the order of, or smaller than, about 30 km. In order 
to facilitate a mutual geostrophic adjustment between the currents and the density field, 


the robust diagnostic method is introduced in the next section. 
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Figure 20. U Component Velocity at Time 72 Hours and Depth of 65 Meters for 
Density Diagnostic: Contour interval 20 cm/s. Dashed lines are 
negative values. 


C. THE ROBUST DIAGNOSTIC AND ITS EFFECTS 

The robust diagnostic method, as described in Chapter II, allows the currents and 
the o, field to mutually adjust to each other. The 0, field is damped back to its initial 
value on a time scale that is fast compared to advection but slow compared to the time 
needed for geostrophic adjustment. As stated above, this type of damping allows the 


sinaller scale features that have short time scales to mutually adjust, but it holds the larger 
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scale features, that might otherwise advect through a sizeable part of the domain 
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The advective time scale is: 





(15) 
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Figure 22. Resultant Horizontal Velocity Vectors at Time 72 Hours and Depth of 
65 Meters for Density Diagnostic: Vectors are plotted at every 4th 
gridpoint. 


where L and U are the length scale and the velocity scale, respectively. In this case, L 
was chosen as 50 km and U was chosen as 50 cm/s which yields a t, of approximately 
10° seconds or 1 day. 


The geostrophic adjustment time scale is: 


yap (16) 
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Figure 23. W Component Velocity at Time 72 Hours and Depth of 90 Meters for 
Density Diagnostic: Contour interval .2 cm/s. Dashed lines are 
negative values. 


where f is the coriolis parameter. For this case, t, is approximately 10* seconds or 1/10 
of a day. 

The damping time scale, therefore, should be much longer than 1/10 of a day and 
much shorter than 1 day to produce the desired effect. Clearly. it is not possible to 


satisfy both of these constraints. For this study, | day (t,) was chosen as the damping 
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time scale to maximize the time over which geostrophic adjustment can occur. Since the 


damping time is the same as the advective time, and not smaller, some advective effects 
will be unavoidable. 

Time series of Og, and the u, v, and w components of velocity are used to illustrate 
the adjustment process with the robust diagnostic (Figures 24, 25, 26, 27, 28, and 29). 
The time series of 6, for map points (172,96) and (160,92) and a depth of 65 meters 
shown in Figures 24 and 27 indicate a downward trend with time of 6, occurring at both 
points for the robust diagnostic method. 

The time series of the component velocities show an initial adjustment as in the 
density diagnostic case followed by a series of long oscillations (Figures 25, 26, 28, and 
29). These long oscillations seem to begin at about 30 hours for the u and v component 
velocities, and somewhat later for the w velocity component. The oscillation is 
particularly flat for the v velocity component at map point (160,92). The long oscillations 
at both points are probably caused by advective effects at the longer integration times, due 
to using | day (t,) as the damping time scale. In other words, the longer scale synoptic 
features are not being held perfectly constant. At the model grid point shown, there is 
sinking motion (w < 0 in Figure 25) on the longer time scale. The resulting vertical 
advection of the mean stratification causes the density (6,) to decrease with time at 65 
m (Figure 24). The model is losing its strictly diagnostic character. This trade-off is, 
however, desirable in terms of other features of the robust diagnostic as discussed below. 

Comparisons between the u and v velocity component time series for the density 


diagnostic (Figures 13 and 15) and those of the robust diagnostic (Figures 26 and 29) 
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Figure 24. Time Series of o, at Map Point (172,96) and Depth of 65 Meters for 
Robust Diagnostic. 


indicate that the adjustment in horizontal velocity occurs faster for the robust diagnostic 
case. The short oscillations which are characteristic of the density diagnostic method are 


not present in the robust diagnostic case. Whereas the density diagnostic has not 
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Figure 25. Time Series of W Velocity Component at Map Point (172,96) and 
Depth of 90 Meters for Robust Diagnostic. 


completely adjusted by 72 hours, the robust diagnostic has, discounting the long wave 


oscillations which are believed to be advective rather than an initial adjustment. 
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Figure 26. Time Series of U and V Velocity Components at Map Point (172,96) 
and Depth of 65 Meters for Robust Diagnostic. 


Comparisons between the density diagnostic and robust diagnostic time series of the 
w coinponent velocity may be misleading (Figures 12, 14, 25, and 28). The robust 


diagnostic w component velocity time series seem to be oscillating much more than their 
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Figure 27. Time Series of o, at Map Point (160,92) and Depth of 65 Meters for 
Robust Diagnostic. 


density diagnostic counterparts. The robust diagnostic w component velocities are, 
however, an order of magnitude smaller than those of the density diagnostic. Therefore. 


what may seem like larger oscillations are not. The order of magnitude difference 
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Figure 28. Time Series of W Velocity Component at Map Point (160,92) and 
Depth of 90 Meters for Robust Diagnostic 


between the robust diagnostic and the density diagnostic w component velocities is 


discussed further below. 
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Figure 29. Time Series of U and V Velocity Components at Map Point (160,92) 
and Depth of 65 Meters for Robust Diagnostic. 


The 0, field is allowed to change over time in the robust diagnostic. The time 
series in Figures 24 and 27 show this change for two points in the model domain. The 


initial O, field at 65 meters depth for the robust diagnostic case is the same as that for the 
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density diagnostic (Figure 16). The corresponding field at 72 hours is shown in Figure 


30. The decrease in 0, with time shown at the two points in the time series (Figures 24 
and 26) is not uniform throughout the model domain. In general, areas of downwelling 
coincide with areas where 0, decreases between 0 and 72 hours, and areas of upwelling 
coincide with areas where 0, increases between 0 and 72 hours (Figures 31 and 35). 
Downwelling brings warmer, less dense water to the 65 meter level and upwelling carries 
colder, denser water to that level. 

In addition to the larger scale changes in 6, due to vertical advection, there has also 
been a tendency for adjustment on the smaller space scales to produce a smoother 0, field 
at 72 hours. Some of the smoothing effect is accomplished by the biharmonic eddy 
diffusion and some by means of the dynamic adjustment. A scale analysis of the 


thermodynamic equation in the model indicates that the parameter r, 


peace (17) 


FL . 
where A = biharmonic diffusion coefficient, gives the ratio of the two effects (i.e. the 
ratio of the biharmonic diffusion term to the change of 6, on the time scale of adjustment 
(f')). Using the value of A from Table 2, we find that r = 5 for L (wavelength) = 10 km, 
while r = .3 for L = 20 km and r = .06 for L = 30 km. Thus, the smoothing effect is 
primarily accomplished by the dynamic adjustment for disturbance wavelengths equal to 
or greater than about 20 km. Biharmonic eddy diffusion dominates for wavelengths of 


10 km or smaller. 
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The initial u and v component velocity fields, and the resultant velocity vectors at 


65 meters depth for the robust diagnostic, are the same as those for the density diagnostic 
(Figures 17, 18 and 19). The corresponding fields at 72 hours are shown in Figures 32, 
33, and 34. Comparisons between the u and v velocity component fields at 72 hours for 
the density diagnostic and the robust diagnostic show that the robust diagnostic fields are 
smoother (Figures 20, 21, 32 and 33). The smoothness of the u and v component velocity 
fields for the robust diagnostic at 72 hours is consistent with the relative smoothness of 
the 0, field for the robust diagnostic at 72 hours as compared to the initial 0, field. If 
the noise in the 72 hour u and v component velocity fields in the density diagnostic case 
is caused by the effects of small scale noise present in the initial objectively analyzed 
fields, the robust diagnostic adjustment eliminates these effects. Also, note that the 
strongest u and v component velocities at 72 hours for the density diagnostic are 
considerably stronger than those for the robust diagnostic (Figures 20, 21, 32, and 33). 
The resultant vector fields at 65 meters depth shown in Figures 22 and 34 reflect this 
difference. The density diagnostic vector field has a maximum velocity over 100 cm/s 
(Figure 22). The robust diagnostic vector field, on the other hand, has a maximum 
velocity of around 75 cm/s, which is closer to the maximum velocities stated above 
observed by Huyer et al. (1990) at 50 meters depth during CTZR1 (Figure 34). 
Comparing the w component velocity field for the robust diagnostic solution at 72 
hours (the initial w is zero by definition for the robust diagnostic as well as the density 
diagnostic) with the corresponding field from the density diagnostic solution shows a 


dramatic difference. As noted in connection with Figure 28, the maximum upwelling and 
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Figure 30. Sigma-0 Field at Time 72 Hours and Depth of 65 Meters for Robust 
Diagnostic: Contour interval .2 kg/m’. 


downwelling velocities for the robust diagnostic are an order of magnitude smaller than 
those for the density diagnostic (Figures 23 and 35). For the robust diagnostic, the time 
series of w component velocity at the two points (Figures 25 and 28) indicate that 
downwelling is occurring at both points at 72 hours at the rate of approximately .02 cm/s 
(17 m/d). For the density diagnostic, on the other hand, the time series of the w 


component velocity (Figures 12 and 14) indicate that upwelling is occuring at both points 
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Figure 31. Sigma-@ Difference Field (Time 72 Hours - Time 0) at Depth 65 Meters 
for Robust Diagnostic: Contour interval .2 kg/m’. Dashed lines are 
negative values. 


at 72 hours at the rate of 200 to 700 m/d. 
: As indicated above, these vertical velocities for the density diagnostic were much 
. higher than the values observed by Washbum et al. (1990) and Kadko et al. (1990). The 
robust diagnostic vertical velocities as represented by the time series at the two points and 


the w field at 90 meters depth at 72 hours (Figures 25, 28, and 35), on the other hand, 
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Figure 32. U Component Velocity Field at Time 72 Hours and Depth of 65 Meters 


for Robust Diagnostic: Contour interval 20 cm/s. Dashed lines are 
negative values. 


are of the order of 20 meters per day, which are similar to the above stated observed 
values of 6 to 25 meters per day. 

The advantages of the robust diagnostic as compared to the density diagnostic are 
now clear. The robust diagnostic solution comes to a new balance much faster than does 


the density diagnostic solution. In addition, at this new state of balance the robust 
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Figure 33. V Component Velocity Field at Time 72 Hours and Depth of 65 Meters 
for Robust Diagnostic: Contour interval 20 cm/s. Dashed lines are 


negative values. 


diagnostic velocity fields, horizontal and vertical, are much closer to the observed values 
than the density diagnostic fields. For this reason, the balance of the discussion focuses 
on the robust diagnostic version of the model. 

Two minor drawbacks of the robust diagnostic require further investigation. The 


first, as mentioned above, is the long oscillations which occur at the longer integration 
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Figure 34. Resultant Horizontal Velocity Vectors at Time 72 Hours and Depth of 
65 Meters for Robust Diagnostic: Vectors are plotted at every 4th 
gridpoint. 


times of the model. These oscillations could be damped by shortening the damping time 
scale; thus, reducing the likelihood of the long wave synoptic features adjusting with time. 
This shortening of the damping time scale could, however, reduce some of the desirable 
geostrophic adjustment. Sensitivity studies, therefore, are needed to determine the 


optimum damping time scale. 
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Figure 35. W Component Velocity at Time 72 Hours and Depth of 90 Meters for 
the Robust Diagnostic: Contour interval .02 cm/s. Dashed lines 
are negative values. 


The second drawback is the boundary noise seen in the w component velocity field 
at 72 hours (Figure 35). The noise is due to the convergence at the boundaries of inertia- 
gravity wave energy that has rapidly dispersed away from the interior of the domain 
during the geostrophic adjustment process. This noise is not as noticeable in the 


horizontal velocity fields because it is primarily confined to the (smaller) divergent part 
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of the currents. The most effective method of reducing the boundary noise would be to 
increase the strength of the sponge, as described above, or to introduce a divergence 
damper in the (u,v) equations at the boundaries. These drawbacks, however, do not 


hamper the results of the interior solution as presented here. 


D. INTERPRETATION OF THE VERTICAL VELOCITY FIELD IN THE 
ROBUST DIAGNOSTIC SOLUTION 


One of the most significant results of the model generated currents is the predicted 
field of vertical velocity. The vertical velocity is important in determining the fate of 
heat, salt and biological material that is advected offshore in squirts and jets in the coastal 
region. Comparing the w field in Figure 35 with the current field in Figure 34 indicates 
that the offshore jet is associated with sinking motion, while the onshore part of the jet 
is associated with rising motion. In this section an attempt is made to understand this 
result in terms of a quasi-geostrophic vorticity budget. 


A simplified form of the vorticity equation is 


a ee pv: 18 
moe AVY, (18) 


where Bis is the total derivative of the relative vorticity and -(¢+£)V- 7 is the 


divergence term (Holton, 1979; Pond and Pickard, 1983). The variation of the coriolis 
parameter has been neglected for this scale of motion. 
Amone et al. (1990) and Bower (1991) describe the divergence/upwelling and 


convergence/downwelling phenomena as observed in the Gulf Stream system. The 
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schematic shown in Figure 36 briefly summarizes the interpretation by these authors. The 
above equation aids in analyzing the schematic. 

As seen from equation {18), if divergence (V- V>0) is present, then the vorticity 
tendency of the fluid particle is negative, implying a tendency for anticyclonic motion. 
If convergence (V- V<0) is present, then the vorticity tendency is positive, implying a 
tendency for cyclonic motion (Figure 36). Conversely, if a fluid particle flows from a 
region of negative vorticity (a ridge) into a region of positive vorticity (a trough), 
experiencing a positive vorticity tendency, then, according to equation {18}, it must flow 
through a region of horizontal convergence. 


From continuity and equation (18), the following relationship is true 


dl _ Ow {19} 


From equation {19} with w = 0 at z = 0, a negative vorticity tendency implies upwelling 


(ow 
Oz 


<0), and a positive vorticity tendency implies downwelling ( OW59) (Pond and 


‘Oz 


Pickard, 1983). In Figure 36, the vertical upwelling centers are shown upstream of a crest 
in an area of divergence and downwelling centers are shown downstream of a crest in an 
area of convergence. 

The resulting vertical velocity in the robust diagnostic model (Figures 35 and 37) 
agree closely with the patterns of Amone et al. (1990) shown in Figure 36. At 72 hours, 
areas of downwelling are located downstream of a ridge in the offshore jet in the south 


central section of the model domain at 65 meters depth. The main area of upwelling ts 
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Figure 36. Convergence and Divergence Patterns in a Typical Meander: Vertical 
upwelling centers are the dots, and downwelling centers are the 
crosses(from Amone et al, 1990). 


located upstream of the ridge where the flow in the jet is directed onshore (Figure 37). 
The vorticity field at the same time and depth (Figure 38) supports the above 
interpretation. The strongest sinking motion in Figure 35 occurs in the region of largest 
vorticity tendency. This is where fluid particles are flowing across the vorticity "front" 
south of the offshore jet. 

The initial pressure and vorticity fields are shown to illustrate the evolution of these 
fields with time (Figures 39 and 40). Note the smoothing which occurs in the vorticity 
field, in particular, between the initial and 72 hour times (Figures 38 and 40). This 
smoothing is the result of the geostrophic adjustments caused by the robust diagnostic as 


discussed above. However, the basic pattern of positive vorticity to the left of the jet 
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(looking downstream) and negative vorticity to the right of the jet is similar to the 
vorticity pattem of 72 hours. 

Equation {19} can be used to estimate the magnitude of upwelling and downwelling 
expected in the CTZ through scale analysis. These arguments also were used by De Jesus 


(1990). The right hand side of equation {19} can be scale analyzed as 
OW pln {20} 
I f H 


where w,,, is the vertical velocity at depth H, and H is the depth of interest (65 meters). 


The left hand side of equation {19} becomes 
Z 
qu" {21} 


when scale analyzed, where U and L are the velocity scale and length scale, respectively. 


Equating the right hand and left hand sides, and solving for w gives 
wR (22) 


U 


Fr is the Rossby number. 


where R,= 


Using the scaling estimates from the advective and geostrophic time scale arguments 


described above 
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Figure 37. Pressure Field at Time 72 Hours and Depth of 65 Meters for Robust 
Diagnostic: Vertical velocity pattem from Figure 35 shown with, 


+ = upwelling, - = downwelling. 
and 
f ~ 10% sec. 
Therefore, 
R - 0.1 


60 





Cruise CTZR1, O6—12 July 1988 


VORT hour 72, depth65 m 


Ye 
"WG " 


= 
E 
x 
o 
(3) 
Cc 
je} 
wo 
2 
a 


. 


. 
Se od 


[ astaaotade 


168 32 48 64 80 96 
Distance (km) 





Figure 38. Vorticity Field at Time 72 Hours and Depth of 65 Meters for Robust 
Diagnostic: Contour interval 1 X 10° sec’. Dashed lines are 
negative values. 


and 


w ~ Ol cm sec"! 


or about 6 m day", which is similar to the lower observed value of Washbum et al (1990) 


and Kadko et al. (1990) and similar to the result shown in Figure 35. These scaling 
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Pressure Field at Time 0 and Depth of 65 Meters for Robust 
Diagnostic: Contour interval 2.5 cm dynamic height. Dashed lines 
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Figure 39. 


arguments indicate that the observed values and the robust diagnostic model values are 
probably characteristic vertical velocity values for the CTZ region. The interpretaion of 
the w field, as that due to the convergence/divergence pattem that changes the vorticity 
of fluid particles as they move through a meandering jet with alternating regions of 


positive and negative vorticity, is probably also characteristic of the CTZ region as it is 
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Figure 40. Vorticity Field at Time 0 and Depth of 65 Meters for Robust 


Diagnostic: Contour interval 1 X 10° sec’. Dashed lines are 
negative values. 


of the Gulf Stream (Bower, 1991). 


E. ADCP COMPARISONS 
The resultant velocity field for the robust diagnostic at tume 72 hours and 65 meters 
depth shows a relatively smooth southward to southwestward flowing jet through the 


center of the model domain (Figure 34). The magnitude of the velocity following the jet 
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is fairly steady, and it’s direction changes smoothly. The ADCP velocity vectors for 


CTZRI at 50 meters depth (Figure 41), on the other hand, indicate the instantaneous 
observed velocity field is not quite so smooth as indicated in the model. The general 
shape of the jet in the observed field is reflected in the model field. The magnitude of 
the velocities in the observed field, however, change spatially to a much greater degree 
and with more irregularity than in the model. This is probably to be expected since the 
ADCP data represent instantaneous currents, while the model currents represent 
geostrophically balanced currents estimated from a fairly coarse sample grid. Also, the 
maximum observed velocities are in the southwestern comer of the cruise domain where 
the jet is farthest offshore. In the model, maximum velocities for the jet are indicated in 
the northeastern part of the domain, closer to the shore. 

The smoothing which occurs because of the robust diagnostic seems to eliminate 
some of the physical small scale disturbances in the velocity field as well as those 
introduced by the objective analysis as discussed above. Vertical cross-sections of the 
velocity field (below) also indicate that the model does not represent well the smaller 
scale characteristics of the observed velocity field. This is to be expected, to some extent, 
because the model currents represent a balanced state, while the ADCP data are 
instantaneous currents. 

Vertical cross-sections of normal velocity from the model domain are compared to 
the corresponding ADCP normal velocity cross-sections from the cruise (Figure 42). The 
cross-sections corresponding to cruise leg B in the northeastern comer of domain show 


that the observed jet (Figure 43) is smaller in vertical and horizontal extent, and has a 
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Figure 41. The ADCP Velocity Vector Field at 50 Meters for CTZRL. 


lower velocity than the jet developed by the model (Figure 44). The position of the jet 


in the model is also shifted to the southeast. The observed position of the jet in cruise 
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leg C (Figure 45) corresponds closely with that developed by the model (Figure 46), but 


the observed intensity of the jet is not as great as that developed by the model. The 
cross-sections for cruise leg D have similar characteristics to those of cruise leg C(Figures 
47 and 48). Cruise legs E, F and G cross-sections, however, show an increasing 
agreement in the magnitude as well as the position of the observed jet and the modeled 
jet (Figures 49, 50, 51, 52, 53 and 54). In general, note that the model velocity cross- 
sections for all legs are much smoother than those for the ADCP velocities. Smaller scale 
disturbances do not appear in the model cross-sections because of the balanced nature of 
the currents that are represented, and the lower resolution of CTD stations used to 
initialize the model. The discrepancy between model and ADCP data below the surface 
jet (i.e. below 200 m), however, cannot be attributed to smaller scale imbalanced flows 
in ADCP data because the currents at that depth should be in dynamic balance. 

There are several possible reasons why the ADCP and model currents are in closer 
agreement farther offshore. Topographic effects, which are neglected in the model, could 
produce the increasing agreement between observed velocities and modeled velocities 
with progression from northeast to southwest through the model. Boundary effects in the 
model could also play a part, especially in the northeast comer of the model domain 
where inflow of the jet is occurring. In addition, the boundary probably affects the model 
cross-section corresponding to cruise leg G (Figure 54). Note the reverse flow which 
occurs in the extreme southeastern part of this leg which does not occur in the observed 
velocity field (Figure 53). More rigorous statistical studies, and robust diagnostic studies 


with topography, wind and heat flux forcing, and barotropic currents included, are 
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Figure 42. Model Domain Cross-Sections, Cruise CTZR1: Solid lines indicate 
model cross-sections. Diagonal rows of dots are cruise legs A 


through G from northeast to southwest. 
required before a better interpretation of the comparison between ADCP data and model 


currents can be made. 
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Figure 44. Model Velocity Cross-section Cruise CTZR1 Leg B: Dashed lines are 
northeastward flow. Northwest is to left. Contour interval 5 cm/s. 
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Figure 45. ADCP Velocity Cross-section Cruise CTZR1 Leg C: Conventions and 
contour interval are same as in Figure 43. 
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Figure 46. Model Velocity Cross-section Cruise CTZRI Leg C: Dashed lines are 
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northeastward flow. Northwest is to left. 
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Figure 47. ADCP Velocity Cross-section Cruise CTZR1i Leg D: Conventions and 
contour interval are the same as in Figure 43. 


12 





oO 
oO 
22) 
2 
=) 
s 
N 
ice) 
© 
a 
N 
i 
O 
@ 
Y 
e 
O 


norm field at hour 72 from (112,156) to (188,80) 





72 80 88 36 104 


32 40 48 56 & 
NW-SE Distance (km) 
Figure 48. Model Velocity Cross-section Cruise CTZRi Leg D: Dashed lines are 


16 24 


0 





Contour interval 5 cm/s. 


northeastward flow. Northwest is to left. 





73 


a) 
SS 
wn 
© 
~~ 
> 
Es 
O 
(@) 
4 
(x) 
> 
OQ. 
O 
A 
<< 


CTZ88/GRID IE 


7704130 


= \o° 


== 


008 
(qp) aanssatd 


oor 


0 


0 


G 





Figure 49. ADCP Velocity Cross-section Cruise CTZR1 Leg E: Conventions and 


are contour interval are the same as in Figure 43. 
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Figure 50. Model Velocity Cross-section Cruise CTZRI1 Leg E: Dashed lines are 
northeastward flow. Northwest is to left. Contour interval 5 cm/s 
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Figure 51. ADCP Velocity Cross-sectiun Cruise CTZR1 Leg F: Conventions and 
contour interval are the same as in Figure 43. 
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norm field at hour 72 from (44,120) to (132,32) 
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Figure 52. Model Velocity Cross-section Cruise CTZR1 Leg F: Dashed lines are 
northeastward flow. Northwest is to left. 
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Figure 53. ADCP Velocity Cross-section Cruise CTZR1 Leg G: Solid lines are 
southwestward flow. Northwest is to left. Contour interval 5 cm/s. 
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Cruise CTZR1, 06-12 July 1988 


norm field at hour 72 from (12,92) to (100,4) 





0 8 16 24 32 40 48 56 64 72 80 88 96 104 12 120 
NW-SE Distance (km) 


CTZRI1 Leg G: Arrows are 


Figure 54. Model Velocity Cross-section Cruise 
boundaries of corresponding ADCP cross- 


Northwest is to left. 


section. Dashed lines are northeastward flow. 


Contour interval 5 cm/s. 
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IV. SUMMARY AND CONCLUSIONS 

The "interface matching" formulation of determining the deep solution for the 
disturbance 0, gives reasonable results when input into the primitive equation model used 
in this study. No comparisons, however, were made between these results and those 
which would be obtained using the "mean conserving” or the "energy conserving” 
formulations. A deep profile comparison of the “interface matching", the "mean 
conserving", and the "energy conserving” formulations with observed data is needed to 
determine which method gives results which most closely resemble the observed vertical 
ocean structure. In addition, a comparison is needed between the results obtained by the 
exact method and those obtained by the linear least squares method of solving the 
minimization problem set up between equations {1} and (2) as described in Chapter II. 

The robust diagnostic version of the model gives better results than the density 
diagnostic version of the model for this case. The robust diagnostic version of the model 
allows a mutual dynamic adjustment to occw between the 6, and velocity fields. 
Consistent with the theory of geostrophic adjustment, the smaller scale features in the 0, 
field are altered the most by this adjustment process. The improvement of the robust 
diagnostic over the density diagnostic versions, therefore, will probably hold true for the 
other CTZ cases as well. The accuracy with which the robust diagnostic version of the 
model diagnoses the magnitude and pattern of the observed vertical velocity associated 


with jets in the CTZ has not been accurately determined. However, the model results that 
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indicate horizontal convergence and sinking motion of the order of 5-20 m/day all along 
the offshore directed part of the jet are certainly consistent with the few observed 
estimates that are presently available. A sensitivity study to determine the optimum 
damping time scale for the robust diagnostic for reducing advective effects is desirable, 
but, as stated in Chapter III, the fact that advective effects were apparent in this particular 
robust diagnostic version did not adversely affect the results. 

Comparisons of ADCP with robust diagnostic model velocity vertical cross-sections 
indicated that small-scale features of the observed velocity field were not reflected in the 
model velocity field. This smoothing of small-scale features by the robust diagnostic was 
seen in all fields, horizontal as well as vertical, The smoothing eliminated some 
undesirable noise resulting from the multi-layer objective analysis scheme, but also 
eliminated some physical features. 

The vertical cross-section comparison between ADCP velocities and robust 
diagnostic model velocities also showed that agreement between the two increased from 
northeast to southwest in the model domain. In addition, ADCP velocities were higher 
in the jet when it was farthest offshore in the southwest corner of the model domain while 
model velocities were higher in the jet where it was closer to shore in the northeast 
comer. A hypothesis for this difference is the model’s neglect of topographic effects, 
wind and heat flux forcing. An interpretation of robust diagnostic model results with 


these effects added is required to test this hypothesis. 
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